Casimir Energy For a Massive Dirac Field in One Spatial Dimension: A Direct 

Approach 
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In this paper we calculate the Casimir energy for a massive fermionic field confined between two 
points in one spatial dimension, with the MIT Bag Model boundary condition. We compute the 
Casimir energy directly by summing over the allowed modes. The method that we use is based on 
the Boyer's method, and there will be no need to resort to any analytic continuation techniques. 
We explicitly show the graph of the Casimir energy as a function of the distance between the points 
and the mass of the fermionic field. We also present a rigorous derivation of the MIT Bag Model 
boundary condition. 
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I. INTRODUCTION 

More than 60 years have passed from the time when H. 
B. G. Casimir published his famous paper [l[ on what has 
been called the Casimir effect ever since. It remained rel- 
atively unknown for over two decades, but from the early 
70s this effect has attracted much attention. The Casimir 
energy is a pure quantum effect with macroscopic man- 
ifestations. Generally, the Casimir energy is defined as 
the difference between the vacuum energies in the pres- 
ence and the absence of any external boundary condi- 
tions or background fields. Both of these energies are in 
general infinite. However the difference between the two 
has almost invariably been calculated to be finite. The 
Casimir effects have been calculated for a variety of fields, 
geometries, number of spatial dimensions, and boundary 
conditions (for a review see @, H, 0, [f| ) . 

Recently the Casimir effect has been studied in connec- 
tion with many physical phenomena. For example this 
effect has been studied in the context of the phenomeno- 
logical chiral bag models of the nucleon d, 0] ■ In such 
models the bag constant, B, is an input parameter to 
the theory. As is well known, this constant is added to 
the Lagrangian density in order to balance the outward 
pressure of the quarks by the inward vacuum pressure 
B on the surface of the bag H . This constant can be re- 
lated to the Casimir energy Q . The Casimir effect for the 
String and the Superstring leads to string theories having 
critical dimensions. The string is a finite two-dimensional 
system with an infinite phonon spectrum, the sum of the 
zero-point fluctuations of which leads to exactly the same 
calculation as in the original Casimir effect [9] . More- 
over, the presence of the Casimir effects in many differ- 
ent phenomena in condensed matter and laser physics 



* Present Address 
Semnan, Iran 
Electronic Address: 
m-valuyan@sbu . ac . ir 



Semnan Branch, Islamic Azad University, 
m . valuyanOsemnaniau . ac . ir ; 



have been established both theoretically and experimen- 
tally EE [ll| . For example, the Casimir force is important 
in the development of microtechnologies that routinely 
allow control of separation between bodies smaller than 
1/im EE EE 0|- An interesting application of the field 
theoretical models with compact dimensions recently ap- 
peared in nanophysics. The long-wavelength description 
of the electronic states in graphene can be formulated in 
terms of a Dirac-like theory in three dimensional space- 
time with the Fermi velocity playing the role of speed of 
light 0, EE EE EE EE HI . Single- wailed carbon nan- 
otubes are generated by rolling up a graphene sheet to 
form a cylinder and the background space-time for the 
corresponding Dirac-like theory has topology R 2 x S 1 [2l[ . 

In order to compute the Casimir energy, physically rel- 
evant boundary conditions must be imposed. For exam- 
ple the Casimir energy for an electromagnetic field is usu- 
ally calculated with the boundary conditions appropriate 
for conducting boundaries as in @, HI, 13, IE HE HE HE HI] ■ 
The Casimir energy for scalar fields has been investigated 
with the Dirichlet boundary condition [E HI, H IE f2p, HE 
HE HE 1 the Neumann boundary condition [E U the 
mixed boundary condition [3l[ , and the Robin boundary 
condition [IE Hi H3]. 

From this point on we concentrate on the Dirac field. 
Since the Dirac field obeys a first order differential equa- 
tion of motion, it is impossible to use the aforementioned 
boundary conditions. Moreover, boundary conditions arc 
in general more disruptive to Dirac fields than boson 
fields because the equations of motion are first order [35j j . 
A proper boundary condition for the Dirac field is the 
MIT Bag Model boundary condition. It is usually said 
to imply that there is no flux of fermions through the 
boundary, i.e. if j M denotes the current of the Dirac field 
and n M is the normal unit vector to the boundary, then 
n^j^ = 0. However, it implies an even stronger condi- 
tion which is the absolute confinement of the fermionic 
field. One could confine the fermionic wave function by 
an infinite scalar potential. This model was considered 
by Bogolioubov [36| and later developed as the MIT Bag 
Model by A. Chodos, et al [IE HI] f° r hadrons. It is 



2 



important to mention that in general the computation of 
the Casimir energy for a massive Dirac field turns out to 
be much more difficult than the massless case. There has 
been relatively few studies for the Casimir energy inside 
closed surfaces, and we mentioned a few of them here. 
The Casimir energy has been calculated for a spheri- 
cal geometry for a massive [39| and massless fermionic 
fields [3j, |4|, |40( . It has also been recently calculated for 
a massless fermionic field subject to the MIT Bag Model 
boundary condition confined inside a three dimensional 
rectangular geometry [4l| . 

As usual there are many more studies on the two par- 
allel plates geometry, and we shall concentrate on this 
problem from this point on. The first computation of the 
Casimir energy for the Dirac field was done by Johnson 
in 1975 [iU ■ He computed this energy per unit area for 
a massless fermionic field subject to the MIT Bag Model 
boundary condition between two parallel plates in three 
spatial dimensions. Afterwards, the Casimir energy has 
been calculated for massless fermions in one dimension 
[35| . between two parallel plates using various methods 
in three-dimensional space J3, & 42|, |43[ and in d + 1 di- 
mensional space-time [3 0, [46 | . The first attempt to 
compute the Casimir energy for the massive case in three 
spatial dimensions was done by Mamayev and Trunov 
who managed to find an integral form for this quantity 
and explicitly computed its small and extremely large 
mass limits [47j . The first complete computation of the 
Casimir energy for the massive case in this geometry in 
3 + f and I + I dimensions was done by C. D. Fosco and 
E. L. Losada in 2008 [48J. In their approach, a coupling 
of the bilinear form tptp to a series of regularly spaced 
(5-function potentials with coupling constant g is intro- 
duced, which implements imperfect bag-like boundary 
condition. This method can produce a fermionic propa- 
gator which satisfies the MIT Bag Model boundary con- 
dition when g — 2. However a direct calculation of the 
Casimir energy for this problem has not been presented 
so far, and this will be the subject of this paper. 

In order to calculate the Casimir energy from first 
principles one must sum over the allowed modes. How- 
ever, the vacuum energies in the presence and absence 
of disturbances obtained by these direct sums are in- 
finite. Therefore one has to adopt regularization and 
analytic continuation prescriptions, in order to compute 
the difference between these two divergent quantities. 
When the allowed modes which appear in the summands 
are regular, the major approaches used are: the zeta 
function analytic continuation technique [3, 0, H l49j . 
cut-off regularization [5(1 HH, [H, HH and box subtrac- 
tion scheme along with the Abel-Plana summation for- 
mula [13, [53|. On the other hand, when the allowed 
modes which appear in the summands are irregular, 
the major approaches used are: the contour integration 
method [H| , the Green function formalism |55|. the func- 
tional approach (48| , and the Boyer method |26l [56j | . Un- 
fortunately, most of these techniques do not be lead to 
closed forms for the values of the Casimir energy and 



one has to employ various numerical methods to obtain 
a value for the Casimir energy. 

In this paper we calculate the Casimir energy for a 
massive fermionic field in two dimensional space-time 
with the MIT Bag Model boundary condition, by di- 
rectly summing over the modes. In this problem the 
fermionic modes are irregular and the divergences that 
appear are very severe. Upon using the contour inte- 
gration method, one might encounter some ambiguities, 
mainly due to the severe nature of divergences inherent 
to this problem. Here we use a direct approach which 
does not resort to any analytic continuation techniques 
and is devoid of any ambiguities. Our method is based on 
the Boyer method, which we shall explain in detail [56jj . 
Moreover, in this procedure one can associate a physical 
meaning to the Casimir energy: It is equal to the work 
done in forming the configuration under study from the 
free vacuum. In section II, we solve the Dirac equation in 
one spatial dimension using the MIT Bag Model bound- 
ary condition and we find a transcendental equation for 
the discrete spectrum. The allowed modes for the mas- 
sive case obtained from this equation are not regular. In 
section III we compute the Casimir energy by performing 
a direct sum over all modes of the field. We finally plot 
the values obtained for the Casimir energy as a function 
of the distance between the points for various values of 
the mass. We show that the results for the small mass 
limit converges to the results for the massless case. The 
results obtained in this paper are consistent with those 
obtained in 48]. In Section IV we summarize our results. 
In the Appendix we present a rigorous derivation of the 
MIT Bag Model boundary condition for the most general 
case. 



II. THE SOLUTION OF DIRAC EQUATION 
WITH THE MIT BAG MODEL BOUNDARY 
CONDITION 

In this section we find the solutions to the Dirac equa- 
tion with the MIT Bag boundary condition in one spatial 
dimension. We consider the case where the Dirac field is 
completely free inside the bag: (— | < x < |). We can 
write the spinor ip(x,t) as: 



ip(x, t) = e 



-iEt I f(x) 



SO) 



(1) 



where E denotes the energy eigenvalue of the time- 
independent solution. We choose the following repre- 
sentation for the 7-matrices: 7 = o\, 7 1 = i<73, and 



7 = 77 



02. Then the Dirac equation leads to, 



f{x) = Ce lkx + De 



g(x) = f(Ce tkx + De~ ikx ) + %{Ce lkx - De~ ikx ), 



(2) 



where k = \J E 2 
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FIG. 1: In this figure the real and imaginary parts of the upper and lower components of the lowest two positive energy 
wavefunctions tl>(x) are plotted as a function of x, separately. In parts. (A,B) the real and imaginary parts of the upper and 
lower components of the ground state (ki = 2.0288 and Ei = + m 2 with m = 1) are plotted. The parity of this state is 

positive. In parts. (C,D) the real and imaginary parts of the upper and lower components of the first excited state (&2 = 4.9132 
and E'2 = +v/fc| + m 2 with m = 1) are plotted. The parity of this state is negative. 



The MIT Bag Model boundary condition has been de- 
rived in the AppendixlAl and we show that this bound- 
ary condition can completely confine the spinor fields be- 
tween the boundaries. The MIT Bag Model boundary 
condition in two space-time dimension in our convention 
becomes, 



Applying these conditions we obtain 



3(f) = 0, 



(4) 



= 0. 



(3) 



X— ±T 



Extracting ^ from the condition <?(§) = and insert- 
ing it into /(-^) = 0, we obtain an expression which 
determines the quantized modes: 

kacot(ka) — —ma. (5) 

The above expression is identical to the result obtained 
by Mamayev and Trunov for the component of the mo- 
mentum perpendicular to the plates for a massive Dirac 
field between two parallel plates in three spatial dimen- 
sions [13] . In Ref. [i?} this condition has been obtained 
by using the MIT boundary condition and the property 
of the Dirac spinors that each component also satisfies 
the Klein-Gordon equation. However all of the allowed 
modes along with this equation for the two parallel plate 



geometry in 3+1 dimension can also be easily obtained 
directly, along the same lines as the derivation presented 
in this section. The solutions of the above transcendental 
equation (Eq. ((SJ), denoted by k s , are not regular for the 
massive Dirac fields. Only for the specific case of a mass- 
less Dirac field, the modes are regular {k s a = ( 2s ~ 1 ) 7r ; 
s = {1, 2, 3, • • • }) [1, 3- Since changing the sign of the 
root k s does not lead to a linearly independent solution, 
we consider only the positive roots. 

As expected the MIT bag model boundary condition 
has transformed the spectrum of the free Dirac field, 
which consisted of two continua starting at E — ±m, into 
two sets of discrete states with energies E — ±y/k^ + m 2 , 
where s — {1,2,3,---}. This problem has parity sym- 
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metry and in this representation the parity operation 
is given by Pip(x,t) — <Jiip(—x,t). The energy eigen- 
states automatically turn out to be parity eigenstates. 
The pari ties of t he lowest lying states with energies 
E = ±y/kf + m 2 are ±1 respectively, and the parities 
of the exited states alternate as the absolute value of the 
energy increases. The real and imaginary parts of the 
upper and lower components of the lowest two positive 
energy wavefunctions ip(x) are plotted in Fig. ([T]). Note 
that the values of the wavefunctions on the boundaries 
are non-zero and just outside the boundaries are exactly 
zero. This is another manifestation of the imposition 
of the MIT bag model boundary condition, as explained 
fully in the AppendixlAl 

In order to calculate the Casimir energy for the mas- 
sive case we need to find a relationship between the root 
number s and the wave number k s . For this purpose, 
Eq. (|5|) can be written as: 



for the Casimir energy, 



X s +tan- 1 (^) 



(6) 



where s is the root number, X s = k s a and M — ma. 
Note that the values of the root indices s and the corre- 
sponding wave-numbers k s can be analytically continued 
to any real value. 



III. THE CASIMIR ENERGY 

In this section we calculate the Casimir energy for 
a massive Dirac field between two parallel plates (two 
points) in 1+1 space-time dimensions. In order to obtain 
the Casimir energy, we should subtract the zero point 
energy in the absence from the presence of the bound- 
ary conditions. Therefore, the vacuum energies for both 
cases should be calculated. Since the solutions of the 
Hamiltonian in the absence and presence of the bound- 
aries are complete [H, [59| , the Fermi field operator can 
be expanded in terms of either of these modes, as follows 



E(x) 



00 dk 

— [b k u k (x) + d\v k (x)} 



^[a s n ks (x) + c\v ks (x)] 



(7) 



where we have denoted the positive-energy and negative- 
energy modes in the free case by u k (x) and v k (x), and 
in the case with the boundary condition by /j. ka (x) and 
v ks (x) , respectively. By substituting the expressions for 
the field operator E(z) given in Eq. (O into the gen- 
eral definition of the Hamiltonian operator and using the 
usual anticommutation relations, and evaluating the two 
zero point energies, we obtain the following expression 



E Cas , = <n i JET | n> — <0 | H bce I 0) 

dx ^2(-y/k 2 + m 2 )vl g {x)v ks (x) 

-oo 



s=l 

+°° '•+ 00 adk 



OO 



dx I ^(_ A / fc 2 + m a) u t( a .) 1 , fc ( a .) ) (8) 



OO 



where | 0) and | fl) denote the vacuum states in the 
absence and the presence of the boundary condition, re- 
spectively. We have just shown that we can obtain the 
vacuum energy, and therefore the Casimir energy, by sim- 
ply summing over the negative energy modes without the 
factor of 1/2. This is equivalent to the usual definition 
where one sums over both positive and negative energy 
modes, since our problem possesses particle conjugation 
symmetry along with C, P and T symmetries, separately. 
Integrating over x the Casimir energy becomes, 



^Cas. = -£v^ 2 

s=l 



+oo 



adk 
~2n 



\[k 2 



(9) 



As mentioned earlier both of these vacuum energies are 
infinite. However the difference, which is the Casimir en- 
ergy, is expected to be finite. In many of the techniques 
for calculating the Casimir energy one starts with only 
the expression for the vacuum energy for the problem at 
hand (the first expression in Eq. ([9]) in our case) and re- 
moves the infinities that appear during the calculation 
by using various methods such as analytic continuation 
or simply by hand, this should precisely amount to sub- 
tracting the free vacuum energy which was left out from 
the beginning (the second expression in Eq. ([9]) in our 
case). 

The dependence of the fermionic quantized momenta 
on the mass of the field is one of the distinguishing fea- 
tures of the Fermi field as compared to the bosonic case. 
In Fig. ([2]) we show the wave vectors for two massive and 
a massless fermionic field. However as mentioned ear- 
lier, for m = the wave vectors are evenly spaced and 
are given by k s — [jfzJjl _ Therefore the Casimir energy 
for a massless fermionic field can be easily obtained from 
the first term in Eq. (JSJ) using the zeta function and its 
analytic continuation as follows d, [35[ , 



E, 



(0) 
Cas 



(m = o) 



OO 

E 

s=l 



(s 



Analytic Part 



-C(-i^) 



— 7T 

24a' 



(10) 



As shown in Fig. ([2]) the wave numbers obtained from 
Eq. ([5]) for a massive Dirac field are irregular. In or- 
der to calculate the Casimir energy we use the Boycr 
method [56| , instead of using Eq. © directly, since the 
latter is more prone to ambiguities. These two methods 
for calculating of the Casimir energy are equivalent. Now 
we discuss the Boyer method in detail. In this method we 
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12 : 
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m=0 

m=l 

m=5 



the Fig. ([3]), for example the region al, we calculate the 
vacuum energy as follows, 



sg ) (M) = Um 

A— >0 



1 oo 
n ' * 



(12) 



6 ; 

4 ; 

2 ^ 

: , 



1 



FIG. 2: Fhe plot of the allowed values of the wave-number 
k s , obtained from the roots of Eq. ((5|, as a function of the 
root number s. We have displayed the results for three values 
of the mass m = {0, 1, 5} with a = 1. Note that for m = 
we obtain a straight line, i.e. the roots are equally spaced. 
However, as is apparent from the figure, this no longer true 
for m ^ 0. 



consider two similar configurations: two points with dis- 
tance a and two points with distance b. Then we place 
each of these systems in a box with size L > a,b as 
Fig. ©. Finally we define the Casimir energy as the dif- 
ference between the vacuum energies of these two similar 
configurations as follows, 



of a— ¥oc 



lim [E^-E^] 

L/o— »oo 



(ii) 



where 4°) = + 22?£> and = < T 



2E; 



(0) 



b2 



bl 



b2 



where we have introduced a convergence factor g(Xio s ) — 

e -W x °+ M2 with A — > 0. Now we start the calculation 
for the region al, and the calculations for the other re- 
gions can be done analogously. Since the wave- vectors are 
not regular with respect to the root indices s, in order to 
find an analytical form for the divergence in Eq. (|12[) , we 
employ the simplest form of the Euler-Maclaurin Sum- 
mation Formula (EMSF) to obtain 



4?V) = -ilim 
a A->o 



ds(X 2 + M 2 ) 1/2 g{Xoj s ) 



s=l 



(13) 



° ( 
dsl i 

=1 ^ 



I)^[(X s 2 +M 2 )VVAc s )] 



where [s] is the Floor function. There is a one-to-one 
correspondence between the wave vector k s and the root 
number s, as is manifest in Eq. ([6]). Therefore we can 
change the variable of integration from s to X . Then we 
obtain, 



E$\m) = --lim 
a a->o 



M 



2x1/2 



Sr(Awi) 



dX^X^+M^g^) 

(14) 



X=Xi 



dXls- [s] 



\)^[(X 2 + M^/ 2 g(Xu>)] 



a2 



al 



a2 



FIG. 3: The geometry of the two different configurations 
whose energies are to be compared. The labels al, 61, etc. 
denote the appropriate sections in each configuration sepa- 
rated by points. The upper configuration is denoted by 'B', 
and the lower one by 'A'. 

are the vacuum energies of configurations 'A' and 'B' as 
shown in Fig. ([3]) . Note that upon taking the limits indi- 
cated in Eq. (|lip , one recovers the original definition of 
the Casimir energy given in Eq. ([9]) . From this definition 
one can easily conclude that the Casimir energy is equal 
to the work done on the configuration 'B' to deform it 
to configuration 'A'. For each of the six regions shown in 



Now, by adding and subtracting appropriate terms, we 
can extend the lower limits of all integrals in Eq. (fl4|) to 
zero. We thus have 



x. 



■— lim 

a A-s-o 

ds 



dX^-(X 2 +M 2 )V 2 g(\u>) 
x=a « A 



dX — (X 2 + M 2 ) 1 /2 5(Aw) 

X=Q « A 



1 dX (s-[s]-^±[{X* + M^g{Xu)] 



(15) 



x=o 



The last term in Eq. (fT5|) can be simplified by noting that 
the Floor function [s] = in the indicated domain, and 



6 



integration by parts yields, 
Xl d x( S -[s]- U " 

X=Q 



= ^xi + M^g^ + lMe-™ 



Xl rla 

dX^{X 2 +M 2 yl*g(\u>). 
x=o dX 



A = 1/A. In Fig. (j4]) the values of the Casimir energy 
have been plotted as a function of the distance a for var- 
ious values of m. This plot shows that there is a good 
consistency between the results of the massless case and 
massive ones when m — > 0. This figure also shows the 
rapid decrease in the value of the Casimir energy as a 
function of ma. We should mention that our results are 
in an excellent agreement with the previously reported 
result which was obtained indirectly through the analy- 
sis of the bilinear forms EMl • 



(16) 



Using Eqs. (|15I16|) we obtain 



— lim 

a a->o 



+ £^ dX (s - [s] - I) ± [(X 2 + M*)V*g{\u>)\ 



(17) 

i 

— AM 



Note that only the first term on the right hand side of 
Eq. (flTl) is divergent. Upon substituting the expression 
displayed in Eq. (IT7l) into the definition of the Casimir en- 
ergy given in Eq. dill) , the constant terms automatically 
cancel each other in the limit A — > and by choosing 
appropriate cutoffs on the upper limits of the integrals 
for each region, the divergent integrals cancel each other, 
all due to the box subtraction scheme. Therefore only 
the convergent integral terms remain. The contributions 
of integrals in regions a2, 61 and 62 go to the zero as 
L/b — > oo and b/a — > oo. Therefore our final expression 
for the Casimir energy is, 



Ecas. = Jim 

A->0 



dX (i 



x^[(X 2 + M 2 ) 1/2 3(Ac)] 



(18) 



where s is obtained from Eq. ([6]) . It seems that this 
expression does not have a closed form solution and 
it should be solved numerically. As is apparent from 
Eq. (|18l) , the integrand has an infinite number of discon- 
tinuities due to the presence of the Floor function. First 
the precise positions of the jumps in the integrand have to 
be determined. These jumps precisely correspond to the 
roots of Eq. ([3]), giving the values of the wave- numbers. 
The integrations are done separately for all parts and 
then all of the results are summed. The integration is 
over the continuous version of the wave number, which 
extends to infinity. In order to accomplish this numer- 
ically, we compute this integral up to a cutoff A which 
should eventually go to infinity. Meanwhile, we also have 
to take the limit A — > as indicated in Eq. (|18p . We have 
determined that an optimization occurs precisely when 




FIG. 4: The values of the Casimir energy for the massive 
and massless Dirac fields with the MIT Bag Model boundary 
condition in one spatial dimension are plotted as a function 
of the distance between the points (a). In this figure we have 
shown a sequence of plots for m = {1, 0.1, 0.01, 10 -6 , 0}. It is 
apparent that the sequence of the plots for the massive cases 
converges rapidly to the massless case as m decreases. 



IV. CONCLUSION 

In this paper we have computed the Casimir energy for 
a massive Dirac field with the MIT Bag Model boundary 
condition in one spatial dimension. We have used the 
direct mode summation method in order to compute the 
Casimir energy for this field. For the massless case the 
modes are regular and we use the zeta function analytic 
continuation technique. However, for the massive case 
the modes are irregular and we use the Boyer's subtrac- 
tion scheme. Its worth mentioning that in this technique 
all of the infinities are canceled automatically and there 
is no need to use any analytic continuation. We have 
shown that the massless limit of the massive case pre- 
cisely corresponds to the massless case. Our result for 
the values of the Casimir energy has been obtained nu- 
merically, similar to the previously reported results [48| . 
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Appendix A: A Derivation for the MIT Bag Model 
Boundary condition 

In this appendix we present a rigorous derivation of 
the MIT Bag Model boundary condition for the Dirac 
field inside an arbitrary closed surface S. This boundary 
condition ensures the complete confinement of the eigen- 
states of the Dirac Hamiltonian inside an enclosed area. 
We show that this boundary condition can be obtained 
by coupling the Dirac field to a scalar potential V and 
taking the limit as V — > oo. However we first present 
the reasons why we cannot confine fermions inside an 
enclosed area by the time component of a four- vector po- 
tential Vq. 

Solving the Dirac equation we obtain p 2 = (E — Vo + 
m)(E—Vo — m) which is positive for \E— Vq j > m. There- 
fore we have oscillatory solutions out of the barrier, i.e. 
we have currents of particles and antiparticles. This con- 
tradicts the assumption of complete confinement of the 
fermionic field. This is the well-known Klein's paradox 
57) . Next we use the scalar potential V(x) which, as we 
shall see, does not have this problem. The Dirac equation 
with the scalar potential is: 



[ij^-(m + V(S))]i>(x,t)=0. 



(Al) 



Decomposing the spatial components of 7 M <9 M at the sur- 
face into tangential (t) and normal (n) parts, we have: 



= 7 °9o + 7% + i n d n . 



(A2) 



We choose the potential V(x) to be infinite outside of 
the enclosed area and to vanishes inside. First we take the 
integral of the Dirac equation with the scalar potential 
(Eq. (jAlj) ) from a — eh to a (a small interval inside the 



barrier) where a specifies a random point on the surface 
and n is the normal unit vector to the surface at point a 



hdn{i^d^ — m)ip(x,t) = 0. 



(A3) 



When e — > 0, all of the terms vanish except the one con- 
taining the term 7™<9„. Then: 

i(h ■ 7)(i/)(a) — %p(a — eh)) = =>■ ij){a) = %p(a — en)(A4) 

Second, we change the integration domain. This time we 
integrate from a point just inside the volume (a — en) to 
a point outside of the volume (a + eh). 



hdnii^df, - (m + V(x)))ip{x, t) = 0. (A5) 



This time we have to deal with the infinite potential out- 
side the bag and therefore we cannot neglect the term 
eV . Now using the fact that the Dirac equation demands 
ip(a + eh) = and Eq. (|A4[) one obtains: 

- i(h ■ i)<ip(a) = eVip(a). (A6) 
Multiplying Eq. (|A6Jl from left by (n • 7) yields eV = ±1. 
Since e > and 7 = +00 wc conclude eV = 1. 

It is interesting to note that using the time component 
of a four- vector potential for the confinement purpose we 
obtain an inconsistent result: eV = ±4. 

Inserting eV — 1 in Eq. (|A6|) we obtain the MIT Bag 
Model boundary condition: 



(1 + i(h ■ j))ip(x) 



0. 



Boundary 



(A7) 
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